Theory of tunnelling into a multi-band superconductor: decoherence and interference 
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By an exact formulation of tunnelling into a multi-band superconductor in terms of Green's 
functions, we demonstrate that the multi-band feature of the iron-based superconductors can lead to 
novel interference between Andreev reflections and decoherence effect of quasi-particles to Andreev 
reflections in a tunnelling junction. These effects depend on the relative sign of the gaps for s-wave 
pairing, and the sign of the Fermi velocities also matters for tunnelling along the nodal direction of 
d-wave gaps. Experimental probe of such effects could determine both the pairing symmetry and 
the relative sign of the gaps on different bands. 



Introduction: The discovery of iron-arsenide high-Tc 
superconductors [1] raises a new tide in the research of 
superconductivity. [2] The immediate theoretical question 
to ask is what is the pairing symmetry, for which a group- 
theoretical analysis has been a first guide. [3] As for the 
pairing mechanism, a plausible scenario is that the strong 
spin fluctuations at wave vectors connecting the hole- 
like Fermi pockets near the T-point and the electron-like 
pockets near the M-pockets centered at the corner of the 
reduced Brillouine zone (BZ) can lead to repulsive pair 
scattering interactions between the two sets of pockets, 
and therefore trigger a minus relative sign between gaps 
on the two sets of pockets. [4] The relative sign between 
the gaps is therefore a crucial factor to verify/falsify such 
scenarios. As for intra-pocket pairing symmetry, the sit- 
uation is in hot debate. Angle-resoled photo-emission 
(ARPES) measurements find a more or less isotropic gap 
on each pocket, lending support to s-wave pairing. [5] 
Nuclear magnetic resonance (NMR)[6] and penetration 
depth measurements [7] see indications of nodes in the 
gap function, implying d-wave pairing. This is also sup- 
ported by some specific heat measurements. [8] However, 
more recent NMR and penetration depth experiments 
seem to favor node-less pairing. [9] In any case the relative 
sign is unknown. Tunnelling experiments provide signa- 
tures of both s-wave and d-wave pairing. [10] This may 
be attributed to the sample quality, but we shall demon- 
strate that the interpretation of the tunnelling data re- 
quires extra care due to the underlying multi-band fea- 
ture. The crucial idea is that tunnelling into a multi-band 
superconductor causes interference between the bands, 
an analogue of the well-known effect in multi-lead meso- 
scopic systems. Such an effect may encode the relative 
sign between the gaps, but is absent in, or only partially 
captured by existing theories of tunnelling. [11] Given the 
complex experimental situation, a theory that captures 
the multi-band effect exactly in tunnelling is emergently 
called for. 

In the following we first sketch an exact formulation of 
tunnelling into a multi-band superconductor in terms of 
boundary Green's functions on the leads, which we de- 
rive for practical cases. We then discuss the tunnelling 
behavior in junctions involving iron-based superconduc- 
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FIG. 1: Schematic plot of the junction involving a supercon- 
ducting L-lead, a scattering region and a normal R-lead. Here 
n is the tunnelling direction, x and y are the basis vectors in 
the BZ, and 9 is the angle between a Fermi vector and n. 
The four pieces of the M-pocket at the BZ corners can be 
recombined as a closed circle centered at one of the corners. 



tors with all of the four Fermi pockets. We demonstrate 
that the multi-bands lead to novel destructive interfer- 
ence between Andreev reflections and decoherence effect 
of quasi-particles to Andreev reflections (both involv- 
ing different bands). These effects depend on the rela- 
tive sign of the gaps for s-wave pairing, and the sign of 
the Fermi velocities also matters for tunnelling along the 
nodal direction of d-wave gaps. Experimental probe of 
such effects could therefore determine both the pairing 
symmetry and the relative sign of the gaps on different 
bands. 

Exact formulation: Suppose we are dealing with an N- 
S junction with a superconductor carrying multi-bands 
on the left-hand lead (L), a normal metal on the right- 
hand lead (R), and a scattering center that couples to the 
cv-band of the leads with amplitudes t a . The junction ge- 
ometry is illustrated in Fig.l. In our case a = r^a, 
on the L-lead and a = N on the R-lead. The tunnelling 
current at bias voltage V can be written as, [12] 

I=~( dojTiG r c r L G a c r R a 3 [f(uja Q ) - f(oja - eVa 3 )] 
J duGl^Gl^lfiu + eV) - /(w - eV)\l) 
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where Gc^ is the retarded (advanced) Green's function 



on the scattering center, r r ' 



-i(£ m ' a - S m - r ) with 
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s ™ = EJC) 2 ^™ ^ (for m = L,R and applied sep- 
arately for the retarded and advanced cases) as the self 
energy correction to the scattering center due to cou- 
pling to the leads so that Gj 1 = z<r — ^2 m S m . Here cr's 
are Pauli matrices in the Nambu space, and the doubly- 
numbered subscripts in Eq.(l) refer to the matrix ele- 
ments. Moreover, G a is the Green's function in the a- 
band at the open boundary before attaching the scat- 
tering center. The first term in Eq.(l) describes quasi- 
particle contribution, while the second is due to Andreev 
reflection which is nonzero if T, L contains off-diagonal 
components. Throughout this work z = u) ± i0 + for re- 
tarded/advanced Green's functions, which is left implicit 
for notational convenience. 

Eq.(l) can be reached by intuitive arguments along the 
line of Blonder-Tinkham-Kapwijk (BTK).[13] The ma- 
trix T = G^T L G^T R describes the transmission prob- 
ability [12] for particles/holes from R- to L-lead. For 
electron transmission, the probability flux is T\\v, with 
v the group velocity. Multiplied by the number of modes 
dp/2-Kh within a momentum interval dp, the total flux is 
Tuvdp/h = Tudui/h within an energy interval duo. In 
equilibrium this must balance the counter flux from L 
to R. Under a bias voltage, an electric current arises as 
doj(e/h)Tn[f(uj) — f(uj — eV)\. Similarly, the hole trans- 
mission contributes a current doj(— e / h)T 2 2[f {oj) — f(u) + 
eV)], where — e accounts for the positive charge. The off- 
diagonal parts of T describes incoming electrons/holes 
from R transmitted as outgoing holes /electrons on L. 
This does not contribute to current, and in fact cancels 
out due to charge conservation. On the other hand, one 
can take the Nambu components as generalized channels, 
by which pairing is nothing but channel hybridization. 
On the R-lead, such inter-channel transmission is exactly 
the Andreev reflection. In the spirit of T, the proba- 
bility reads A = G r c A2 ^22^c ^l^fi- This contributes an 
electric current dui(2e/h)Af(ui — eV), where the factor 
of 2 accounts for the charge 2e transported to L during 
each Andreev reflection. Similarly, for incoming holes 
reflected as outgoing electrons, the electric current is 
doj(—2e/h)Af(uj + eV). Collecting everything and in- 
tegrating over lo we arrive at Eq.(l). 

We emphasize that at this stage the formulation is ex- 
act, and all possible interference effects in our case are ex- 
actly reflected in the self-energies that combines the con- 
tributions from all bands coherently. To apply the theory 
we need the boundary Green's function, the derivation of 
which is sketched as follows. 

Boundary Green's functions: Consider an arbitrary 
band in a semi-infinite system with a boundary nor- 
mal n = (cos#„, sin#„) relative to the x direction of 
the BZ (the iron- arsenide-iron direction), as illustrated 
in Fig.l. We assume circular Fermi surface, and near the 
Fermi momenta we fix the conserved transverse compo- 
nent K\\ = Kf s'md, where 9 is the angle between K/ 
and n. For the M-pockets, Q\\, the transverse projection 



of Q = (it, 7r), should be added. We then map the longi- 
tudinal degrees of freedom (near the Fermi angles 9 and 
7T — 9) to a semi-infinite one-dimensional system described 
by the Andreev equations, [13] demanding that the wave 
function vanishes a unit distance (the lattice constant, 
for example) beyond the boundary. These eigen func- 
tions and energies are then used to construct the Green's 
function in a standard manner. Another approach is the 
T-matrix theory for an infinite plane with an infinite bar- 
rier along a line. [14] We get the same results using both 
approaches under the quasi-classical approximation. 

For the s-wave case, the result is independent of 
n: G(9) ~ -AfJ(6)^(A) where Af = K f /2ir\v f \ is 
the normal state density of states (DOS) with vj the 
Fermi velocity, J(9) = 2tt sin 2 (Kf cos 9)/\Kf cos 9\ is 
a Friedel-like factor due to the open boundary, and 
W (A) = {z ap + Agi)/VA 2 -z 2 , in whic h VA 2 - z 2 = 
VA 2 - lo 2 0(A 2 ~lu 2 ) TisignHV^ 2 - A 2 9(u 2 - A 2 ) (for 
retarded/advanced cases) in practice. Here 9(x) is the 
step function. 

For the d-wave case, we consider only two practi- 
cal cases: If n is along the anti-nodal direction, A# = 

Acos26», and we find G{9) AfJ{9)^(A e ) similar to 

the s-wave case, except for the ^-dependence in the pair- 
ing function. 

In contrast, if n is along the nodal direction, Ag = 
A sin 29, an d we find G{9) ~ NJ{9)<S>{A g ) with $(A) = 
(V A 2 — z 2 o~q — As\gti{v f)o~ 2) / z, where sign(wj) = ±1 
if Vf is outward (inward) to the center of the Fermi 
pocket. The pole at u — is a particle-hole mixing bound 
state due to the nodal boundary. [15] Notice that the off- 
diagonal component in $ depends on both the sign of the 
gap and the electron/hole nature of the Fermi pocket, 
leading to interesting consequence for interference in the 
Andreev reflection to be discussed. Finally the Green's 
function in the normal lead can be obtained by simply 
set A = 0. 

The above results are applied separately for G in each 
band, which is then substituted into Eq.(l) for the con- 
ductance calculation. 

Application to junctions with iron-based superconduc- 
tors: We define "f a = A/" a f 2 for the four bands a — 
Fi.2,Mi.2 on the L-side and for a — N on the normal 
R-side, with K a the Fermi radius of the a-pocket. To 
simplify the discussion we assume K a = Kf = tt/2 for 
all bands on the L-lead, and set 7^ = 7 and 7 Q =ri,M; = 
K7/4 for a range of values of k (k = 1/0 corresponds 
to the transparent /tunnelling limit). We also set 7 as 
the unit of energy. For the gap amplitudes on the four 
pockets we set [Af^mJ = Ai and |Ar 2 ,M 2 | = A 2 , with 
Ax = A max and A 2 = 0.6A max in line with existing ex- 
periments. We assume that the gaps on the T-pockets 
are in-phase, so are those on the M-pockets, but the rel- 
ative sign v — ±1 between the T- and M-pockets is left 
free for discussion. 

We consider wide junctions, for which the conduc- 
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tance ag(uj) = dlg/dV is calculated by Eq.(l) (in the 
zero temperature limit) for each conserved transverse mo- 
mentum K\\ = Kj sinO (in which Qu should be added 
for the M-pockets). The channels unmatched by K» 
should be shutoff. The total conductance is then given 
by cr(uj) = jZ^ 2 <7g(ut)K f cos8d8/27T per unit area. 

We first consider tunnelling along n = (1, l)/^/2. Here 
Q|| = 0, and K\\ can be matched to all Fermi sur- 
faces under our assumption, because of which we ex- 
pect novel effects. For the s-wave case, by using G a 
derived above and under the given assumptions, we get 
£ L (#) = -^WE i= i, 2 (^o - A i( j 1 <L,i)/2 v /Af^. 
a) In the in-phase case v = 1, the <j\ -component in T, L 
is left intact, and we expect Andreev contribution. How- 
ever, quasiparticles at \u>\ > A 2 acts as dissipative source 
for (and therefore reduce) the Andreev reflections associ- 
ated with the larger gap. This is a decoherence effect, and 
is seen as the dips in Fig. 2(a) at energies in between the 
two gap edges, b) In the anti-phase case v = — 1, the off- 
diagonal component cancels out in T, L and consequently 
in G c . According to Eq.(l) the Andreev contribution dis- 
appears. This is exactly the anticipated destructive in- 
terference effect among Cooper pairs transported to the 
L-side during Andreev reflections on the R-side. With- 
inin the smaller gap A 2 there are no quasi-particles either 
so that the conductance is completely zero at any level 
of transparency. This behavior is shown in Fig. 2(b). For 
| a; | > A 2 , quasi-particles from T 2 - and M 2 -bands start to 
contribute. Of course, both effects becomes less signifi- 
cant in the tunnelling limit k — > where quasi-particles 
dominate. 

For the ci^-wave case (1, 1) is the antinodal direction. 
According to the previous derivation of G, the structure 
of S L is similar to that in the s-wave case, except that 
Aj should be replaced by Aj cos 26. The decoherence 
effect for v = 1 in Fig. 3 (a) is less visible because of the 
size variation in the d-wave gap function. However, the 
destructive interference effect for v = — 1 in Fig. 3(b) is 
still obvious at k <~ 1. 

For the eL^^-wave case, (1,1) is the nodal di- 
rection. Here the structure of T, L is very dif- 
ferent, £ L (#) = k 7 J(0) E 4 =i, 2 [VA 4 W 2 -^ + 
/\ i (ff)o 2 b-v.\\l2z, where A,(6>) = A;sin(20). The extra 
minus sign before v comes from the fact that the rela- 
tive sign of the Fermi velocities is minus for the T- and 
M-pockets. Therefore we expect destructive interference 
for v = 1 instead of v = — 1. a) The result for v = — 1 is 
presented in Fig. 3(c). For n = 1 the conductance is simi- 
lar to the corresponding line in Fig. 2(a). With decreaing 
k, a central peak shows up and we checked numerically 
that the Andreev contribution dominates, due to the fact 
that the bound state is a particle-hole mixture, b) For 
v = 1, the off-diagonal component cancels out in S L , 
and the results are shown in Fig. 3(d). For k = 1 the con- 
ductance is of broad V-shape. The normal (i.e., in the 
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FIG. 2: Conductance (normalized by the normal state a n ) 
versus the bias energy co — eV for tunnelling along the (1, In- 
direction. The relative sign between the s-wave gaps is v = ±1 
in (a) and (b) respectively. In each panel k = 1, 0.9, 0.1 
from top to bottom, and the curves are vertically shifted for 
clarity (which approach unity on both sides of ui otherwise). 



(a) antinodal junction, v = 1 (b) antinodal junction, v = -1 
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(c) nodal junction, v = -1 (d) nodal junction, v = 1 
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FIG. 3: The same plot as Fig. 2, but for the d-wave case if the 
tunnelling direction (1, 1) is along the antinodal direction for 
(a) and (b), and the nodal direction for (c) and (d). 



absence of particle-hole mixing off-diagonal component) 
bound-state does not show up here because it has zero 
group velocity. In the tunnelling limit n — > 0, the normal 
bound state eventually shows up as a central peak, since 
in this limit the DOS on the L-lead is measured. In addi- 
tion, no coherence peak occur in (c) and (d) even in the 
tunnelling limit, due to the fact that no singularity arise 
at the gap edge in <f>. The absence of coherence peaks 
was possibly observed. [10] 

For tunnelling along h = (1,0), the two sets of pockets 
match that of the normal lead at different momenta, so 
that the T- and M-pockets can be treated independently 
Since within each set of bands the two gaps are in-phase, 
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the effect of multi-gaps can be traced in the previous 
discussion. 

We add without going into details that in-plane point- 
contact on wide boundary is also helpful. Here one 
simply averages the Green's function over 9 first to 
get the on-site Green's function, which is then substi- 
tuted into Eq.(l). There is no momentum matching 
constrain. We see that the off-diagonal part of G sur- 
vives for the s-wave case and the d-wave anti-nodal tun- 
nelling case, in which both interference and decohcr- 
ence effect can occur. In the case of d-wave nodal tun- 
nelling, the off-diagonal part cancels out. Furthermore, 
even out-of-plane tunnelling is useful. Here one sim- 
ply use the on-site Green's function in an infinite plane: 
G a = -7rAf a {[za + A a (#)<7i]/ \J A a (9) 2 — z 2 )g. For s- 
wave gaps, the off-diagonal part still survives and in- 
terference/decoherence occurs, while it cancels out by 
symmetry for all kinds of d-wave gaps and only quasi- 
particles contribute to the conductance. 

We remark that in the real system, the sizes of the 
four Fermi pockets differ slightly. The coupling between 
the scattering center and each band may also differ. 
These factors may weaken the effects discussed so far, 
but should still be discernable. 

Conclusions: We presented an exact formulation of 
tunnelling into a multi-band superconductor, and de- 
rived the relevant boundary Green's functions. Applying 
to iron-based superconductors we demonstrate that the 
multi-band feature can lead to novel destructive interfer- 
ence between Andreev reflections and decoherence effect 
of quasi-particles to Andreev reflections (both involving 
different bands) in a tunnelling junction. These effects 
depend on the relative sign of the gaps for s-wave pair- 
ing, and the sign of the Fermi velocities also matters for 
tunnelling along the nodal direction of d-wave gaps. We 
propose that in-plane wide junction and point-contact 
junction, as well as out-of-plane contact junction can be 
used to probe such effects, and therefore provide decisive 
information on the pairing symmetry and the relative 
sign of the gaps on the various bands. 

While finalizing this work, we become aware of some 
recent interesting theories. [16] 
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tion No. 91009, the Ministry of Science and Technology 
of China (under the Grant No. 2006CB921802 and 
2006CB601002) and the 111 Project (under the Grant 
No. B07026). 
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